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ABSTRACT 



We show how to use 21-cm emission and absorption studies to estimate the heat inputs 
to the neutral gas in low pressure environments, such as in outer disks of spiral galaxies, 
in galactic halos or in intergalactic space. For a range of model parameters we calculate 
the gas kinetic and spin temperatures (Tk and T$) and the relation between Tg and the 
heat input to the gas. We outline the conditions for a "two phase medium" to exist. We 
find that although Tg can be much smaller than Tk, Ts is always > 3 K for column 
densities greater that 5 x 10 18 cm -2 . This excludes the possibility that relevant HI masses 
at the periphery of galaxies are invisible at 21-cm in emission. Therefore sharp HI edges, 
observed in outer disks near column densities Af;~2x 10 19 cm -2 , cannot be "fictitious" 
edges due to a sudden decrease of the 21-cm brightness. The outermost interstellar gas in 
a disk galaxy is more directly affected by external processes and in this paper we estimate 
the intensity of the extragalactic background at energies close to 0.1 keV by comparing 
our theoretical results with HI emission/absorption studies. We take into account the 
possibility that some energy produced in the inner regions affects the energy balance in 
outer regions. We find that in the absence of any other local heat source QSO dominated 
background models are still compatible with the spin temperature limits derived for the 
two best documented HI emission/absorption studies in outer regions. However, if future 
observations should establish that the spin temperature is as high as 1000 K, then relevant 
energy inputs from local sources become necessary. 

Subject headings : Cosmic Background Radiation - Galaxies:ISM - Radio Sources: 21 cm 
Radiation. 



1. INTRODUCTION 



Neutral hydrogen exists outside the stellar region of spiral galaxies as extended disks, 
coplanar or tilted with respect to the innermost luminous disk. Isolated concentrations, 
clouds or plumes have sometimes been detected well outside the HI disk (see for example 
Giovanelli & Haynes 1988 for a review). These blobs may indicate a recent disturbance 
in the HI distribution due to tidal interactions with a companion or to a strong burst of 
activity in the inner disk. In our Galaxy smaller blobs of HI, closer to the luminous disk, 
and known as high velocity clouds (hereafter HVC), have also been detected (Giovanelli 
1980). 

The physical conditions of the gas in the outer regions of spiral galaxies are not constrained 
as they are in the interstellar medium. In the inner disk of our Galaxy, where star formation 
takes place, the following conditions hold: (%)neutral hydrogen column densities are of order 
N HI ~ 2 x 10 20 — 10 21 cm~ 2 , and extragalactic UV or soft X-ray photons with energies 
between 13.6 and 200 eV hardly penetrate this layer; (ii)even without magnetic pressure 
and bulk motion ram pressure, the gas thermal pressure is fairly large, P/k ~ 3000 cm ~ 3 
K (k is the Boltzmann constant); (mjneutral hydrogen exists in a cold phase at the kinetic 
temperature T K ~ 80 and in a warm phase at T K ~ 8, 000 K which fills a non trivial part 
of our ISM (> 30% Brinks 1990; Knapp 1990). The ionization of the gas in the midplane 
might be due to O and B starlight (Kulkarni & Heiles 1988), and the mechanical energy 
input to the gas from O - star winds and from blastwaves produced by supernovae (McKee 
& Ostriker 1977). However these hot stars cannot easily account for the diffuse warm 
component observed at a few hundreds parsecs above the galactic midplane (Reynolds 
1990). 

In outer disks, or in high velocity clouds, stars are absent and the ionizing flux as well 
as the mechanical energy input, from the stellar disk of the galaxy might be very small. 
If chains of supernovae feed galactic fountains which reach the halo directly (Heiles 1987; 
Corbelli & Salpeter 1988; MacLow & McCray 1988; Martin & Bowyer 1990) or if there is 
hydromagnetic wave dissipation (Ferriere, Zweibel & Shull 1988) then the non radiative 
energy input into the outer regions are important. Ionizing flux from extragalactic sources 
is likely to be present due to thermal free-free emission from hot gas, to quasars, or to more 
exotic sources like dark matter decay (Melott, McKay, & Ralston 1988; Sciama 1990) or 
hidden X-UV emitters. The far UV and soft X-ray emission from quasars is of interest 
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here because we mainly consider regions with Njji < 5 x 10 19 cm -2 and photons below 
0.2 keV, which do not affect the inner disk, penetrate these smaller column densities and 
give a large ionization rate. We do not discuss which ionizing fluxes and non radiative 
heat inputs are more likely on physical grounds but we parametrize the overall ionizing 
flux and the non radiative energy input separately and calculate their effects on the HI 
gas. Using HI observations we can then investigate indirectly the energetic environment 
of spiral galaxies. The absence of stars in outer regions makes self-gravity of the gas layer 
an important force which has to be considered, together with dark matter and coronal 
gas, in equilibrium models. Nevertheless we expect a much lower thermal pressure than 
that measured in the inner disk, as suggested by the flaring of the outer disk observed in 
external galaxies and in our own Galaxy (Merrifield 1992). 

HI absorption studies both on outer disks and on HVC, show that only a very small fraction 
of the gas has low enough temperature to be detected in absorption (Corbelli & Schneider 
1990, hereafter CS; Colgan, Salpeter, & Terzian 1990, hereafter CST; Wakker, Vijfschaft, 
& Schwarz 1991; Carilli, van Gorkom, & Stoke 1989; Carilli & and van Gorkom 1992; 
Schneider & Corbelli 1993). However sensitive searches of HI emission at some locations 
where no absorption measurements are possible, indicate that there is a transitional HI 
column density Ni such that appreciable emission is seen for Nb ~ Ni but little emission is 
seen for Nb < iVj. We use the symbol Nb to indicate the column density derived directly 
from the 21-cm brightness temperature; Nb = Nhi when T$ is well above the background 
radiation temperature, Tr, otherwise Nb < Nhi- The transitional column density Ni is 
of the order of 5 x 10 18 cm -2 for high velocity clouds and is of the order of 2 x 10 19 cm -2 
for outer disks (CST; Corbelli, Schneider, & Salpeter 1989) or slightly smaller (Hoffman 
et al. 1993). In this paper we discuss (and then eliminate) the possibility that, due to 
low pressure conditions and to low photoionization rate in outer regions, when Nhi gets 
below N[, the spin temperature strongly deviates from the gas kinetic temperature and, 
while Ts approaches Tr, Nb is depressed with respect to Nhi- Fictitious HI edges would 
then appear if one assumes that Nb — Nhi- The lack of HI absorption lines in outer 
regions already suggests that "subthermal effects" cannot hide appreciable HI masses, but 
we shall show more quantitatively that subthermal effects are not responsible for HI edge 
appearances. A totally different possibility, and a likely explanation for HI edges, is that 
they are sharp HI-HII transitions due to some ionizing flux. This will be analyzed in detail 
by Corbelli and Salpeter (1993, hereafter Paper II). Nevertheless due to the uncertainties 
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in the evaluation of HI masses in extragalactic objects where no absorption measures are 
possible (Giovanelli & Haynes 1991; Schneider et al. 1989) and also to the importance of 
subthermal effects for HI absorption lines, we shall present here a full calculation of the 
spin temperature for several heating inputs outside the galaxy. 

In Section 2 we introduce some models for the extragalactic soft X-ray background ra- 
diation and in Section 3 we display the equilibrium equations and parameters used for 
modeling the gas distribution outside the stellar disk. Section 4 shows results relative to 
subthermal effects in the 21-cm emission line and in Section 5 we derive in more generality 
the relation between the kinetic and spin temperature as a function of the background 
radiation intensity, and of additional non ionizing local heat sources. We outline the con- 
ditions of the HI in a low density medium pervaded with decaying neutrinos in Section 6. 
The last Section summarizes constraints on the energetic environment of spiral galaxies 
which can be set by comparing our theoretical results with 21 cm observations. 



2. LOCAL ENERGY INPUTS AND THE EXTRAGALACTIC SOFT X-RAY 



BACKGROUND 



As discussed in the introduction, there may be energy inputs which are non radiative and 
non ionizing and which may vary from one galaxy to another or from one region to another 
(tidal interaction, hydromagnetic wave heating, fountains, etc.). We shall characterize 
these inputs using a single number, E nr , the heat per hydrogen atom per second due to 
non photoionizing sources in units of eV s _1 , which should be added to the heat input 
from ionizing radiation, E r . The ionization- recombination equilibrium depends on the 
photon spectrum N(E). Below we summarize the known upper or lower limits to the 
extragalactic background radiation at zero redshift and then we will parameterize various 
model assumptions. 

The UV and soft X-ray background has been reviewed by Bowyer (1991) and by McCam- 
mon, & Sanders (1990). Above ~ 2 keV a uniform extragalactic X-ray flux of power law 
form is known to exist (Schwartz 1978): 

^-L = 7.7£ fc -\> 4 ph cm" 2 s" 1 sr" 1 keV" 1 (1) 
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Between 2 keV and 0.5 keV the total X-ray background is about 3x 10 -8 erg cm" 2 s _1 sr -1 
(Shanks et al. 1991; Wu et al. 1991), higher than an extrapolation of eq. (1) to this energy 
range, but the extragalactic contribution is uncertain because of local emission, and it is 
not measurable directly at all below 0.3 keV because of absorption by the local HI disk. The 
space density and luminosity function for quasars are known for redshifts z < 2 (see e.g. 
Hartwick and Schade 1990) and their contribution to the flux between 2 keV and 0.5 keV 
has been estimated. Preliminary analysis of ROSAT deep survey images shows that quasars 
can account for at least 20% of the measured X-ray background below 2 keV (Hasinger, 
Schmidt, & Trumper 1991; Wu, & Anderson 1992) and therefore we consider 5xl0~ 9 erg 
cm -2 s -1 sr -1 as a definite lower limit for the extragalactic background intensity between 
2 kev and 0.5 keV. An extrapolation of eq. (1) down to 1.5 keV yields this minimum 
expected intensity which we consider in examining the possibility of strong subthermal 
effects in emission. 

In this paper we use a few spectral models of the extragalactic background radiation down 
to 0.1 keV. If this background is dominated by quasars, its spectral shape in the absence 
of attenuation is approximately of the form: 

^ = / x E~\j ph cm" 2 s" 1 sr" 1 keV" 1 (2) 
dE keV v ' 

which we shall refer to as the "quasar spectrum". For / ~ 15 the spectrum described by 
eq. (2) corresponds to a flux density ~4x 10 _23 (.E e y/13.6) _1 ' 4 erg cm -2 s _1 sr -1 Hz _1 . 
This flux at zero redshift was suggested by Sargent et al. (1979) and connects the so called 
"blue bump" , observed in QSO spectra, with the cosmic background flux observed above 2 
keV. In a recent paper Madau (1992) estimates a lower value of the integrated UV cosmic 
background from observed QSOs. He takes into account the possible attenuation of the 
flux due to intervening Lyman-a clouds. However above 0.1 keV attenuation is probably 
unimportant and Madau's estimate of the unattenuated flux down to a cutoff energy of 
0.1 keV is given by the quasar spectrum of equation (2) with I ~ 5. 

We shall consider three different cases for flux and energy inputs. For Case A we assume 
the minimum extragalactic X-ray flux in order to compute the most extreme subthermal 
effects. We use a more likely flux for Case B and Case C and show how to use HI spin 
temperature measurements to constrain the extragalactic background flux or the local 
heating just outside disk galaxies. 

Case A : We set E nr =0 and use the flux described by equation (1) with a lower cut- 
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off energy E c = 1.5 keV or E c = 0.5 keV (in this last case about 60% of the observed 
background flux between 2 and 0.5 keV is of extragalactic origin). 

Case B : We consider E nr ^ as an adjustable parameter and use a low quasar flux, 
namely equation (1) with a cut off energy E c = 0.1 keV. 

Case C : We add a soft X-ray component in the energy range E c < E < 1.5 keV, to the 
cosmic spectrum given by eq. (1). The following form describes the spectral law which we 
use in this case : 

{7.7x£"\> 4 E> 1.5 keV 
A E~l v 0.5 < E < 1.5keV 
I E~ e s v E C <E< 0.5 keV 
E <E C 

where A and b are constants determined by requiring dN/dE to be continuous at 0.5 and 
at 1.5 keV and I is an adjustable parameter. In this paper we use E c = 0.1 keV, s = 2.4, 
and E nr = unless stated otherwise. The quasar spectrum in eq. (2) can be represented 
to a good approximation by eq. (3) using s = 2.4 (and / ~ 5 for Madau's intensity). 
In Fig. 1 we show dN/dE for some spectra used. The background spectrum at energies 
close to the Lyman limit will not be of interest in this paper because here we consider 
the heating of slabs with HI column densities £ 0.5 x 10 19 cm -2 . However the complete 
absence of photons below 0.1 keV, is an oversimplification; some photons below 0.1 keV 
will penetrate the gas with low HI column density and increase the ionization rate slightly. 
In Section 6 we discuss Ts and Tk for the case of monochromatic photons with E ~ 14.5 
eV which are produced locally throughout the disk. 



ph cm 2 s 1 sr 1 keV 1 



(3) 



3. CHARACTERISTIC PARAMETERS OF HI GAS OUTSIDE A GALAXY 



The brightness temperature of the 21-cm emission line, Tg, is related to the column density 
of neutral hydrogen along the line of sight, Nhi- 

Tb oc N b = XbN hi , X b = TS ~ TR (4) 

+ s 

where Tr = 2.73 K is the temperature of the microwave background radiation, Ts is the 
spin temperature, defined through the population of the two hyperfine levels of neutral 
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hydrogen (Field 1958). The factor xb is close to unity if the kinetic temperature, Tk is 
well above the cosmic background temperature, and the medium is in LTE (Tg ~ Tk). In 
this case the brightness temperature is directly proportional to Nhi- The spin temperature 
is important for the optical depth at the line center r c which can be determined through 
absorption measurements 

5.14 x 10 19 N HI 

Tc = =- ( 5 ) 

U! Is 

where u is the width of the line at half maximum in km s _1 . Most absorption and emission 
studies give T$ and Njji and not the kinetic temperature of the medium because u might 
not be proportional to \JTk, like if there is turbulence in the medium or the absorption 
line is unresolved. If u is unknown we use the following value: 

uj ~ max(4,0.21 v / 7V) km s" 1 (6) 

based on the line width of the few detected absorption lines (Rubin, Thonnard, & Ford 
1982; Carilli & van Gorkom 1992). Deviations of T$ from Tk are especially important in 
outer regions where low densities slow down the collisional excitation of hyperfine levels 
of H. We calculate Ts for all cases in this paper using equation (15) of Field (1958). We 
consider Lyman-ct pumping through photons generated by H + — e recombinations (Spitzer 
1978), by bound-bound transitions due to thermalized electrons (Spitzer 1978), and by 
bound-bound transitions due to non thermalized electrons coming from H photoionization 
(Shull & Van Steenberg 1985). Following Deguchi & Watson (1985) we set the Lyman- 
a temperature Tl equal to Tk and use the expression of Bonilha el al. (1979) for the 
number of Lyman-o; photon scatters inside the cloud, disregarding the possibility of strong 
velocity gradients. For the collisional de-excitation probability of the hyperfine levels via 
neutral and ionized H atoms, and via electrons, we use the expression given by Smith 
(1966). For T R « T$ « T K the gas is still far from thermodynamic equilibrium but the 
factor xb is very close to unity. This means that subthermal effects show up much more 
easily in absorption rather than in emission measures, because emission depends on the 
ratio between the population of the two hyperfine levels which remains ~ 3. Absorption 
depends on their population difference which can vary by a much larger factor. 
The gas will be considered distributed in a slab, with the vertical extension much smaller 
than the horizontal one. This is a good representation for extended HI disks or flat clouds 
in the outer parts of galaxies. For more spherical blobs our calculation will be indicative 
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but not rigorously correct. The vertical distribution of gas density depends on the ratio of 
dark matter to gas and on any contribution of magnetic fields to pressure (or to buoyancy) . 
These effects will be discussed in detail in Paper II; in this paper we assume that the ideal 
gas law connects pressure, P, to gas kinetic temperature, Tr-. Furthermore we shall use 
here the single slab approximation. For fixed values of the HI column density we compute 
the kinetic temperature and subthermal effects at a height z = z\/2 above the central plane 
such that half the gas mass per unit area lies between —z\/2 and +£1/2 • The attenuation of 
the extragalactic flux due to photoionization of HI, Hel, and Hell, is computed assuming 
that the fractional ionization of Hel and Hell at all levels is the same as at z t / 2 - 
The real pressure will have an additional term due to the external pressure, P ex t, and to the 
compression from a dark matter distribution. Charlton, Salpeter & Hogan (1993) estimate 
that the purely extragalactic value for P ex t/k is less than ~ 4 cm -3 K but some galactic 
coronal gas is likely to increase the effective value of P ex t- For a spherical halo which gives 
a rotational velocity of 100 km/s at R = 15 kpc, the dark matter compression is equivalent 
to about 70 cm -3 K. Since in this paper we use the uniform slab approximation we mimic 
the effects of P ex t, dark matter and eventually of some internal non-thermal pressure in 
one single term, Po- We use the ideal gas law and we consider 0.1 < Po/k < 500 cm -3 K. 
The following expression describes the total pressure we use for a given column density of 
gas (see Spitzer 1942, Ibahez and di Sigalotti 1984): 

where G is the gravitational constant, m# is the hydrogen mass, hu e — 0.1 is the helium 
abundance relative to hydrogen, and N tot is the total (HI+HII) hydrogen column density 
(i.e. twice the gas surface density above the equatorial plane). 

We omit dust grains entirely and we assume that all elements heavier than He have only 
half the solar abundance (unless specified differently) and that there is enough UV light 
below 13.6 eV, from the background or scattered outside the stellar disk, to keep most of 
the carbon, silicon and iron singly ionized. For the range of pressures that we consider this 
is a good assumption (Boyer 1991). Then the heating rate due to photoionization of CI 
does not depend on the UV photon intensity but only on the mean photon energy and we 
use the value Eq = 2eV for the excess energy. In the temperature range that we consider 
(Tr < Tk < 17, 000 K) cooling is dominated by collisional excitation of CII, Sill, Fell by 
electrons and HI atoms at low temperatures, and by excitation of the n=2 level of HI at 
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higher temperatures. However cooling by free-free emission and by HII recombinations are 
considered as well. Po/k > 0.1 ensures the recombination time for HII to be shorter than 
10 10 years. 

For the ionization - recombination equilibrium of HI, Hel, and Hell we use the "on the 
spot approximation" with the photoionization cross section of HI, Hel and Hell as given 
by Spitzer (1978) and Brown (1971). We assume that no photons escape from hydrogen 
recombinations to the n=l level (nor from helium recombinations to n>l) because they are 
absorbed elsewhere by other HI atoms. Also, due to low volume densities in outer regions 
we can assume to a good approximation that all the Hel and Hell recombination photons 
are absorbed by HI atoms (Osterbrock 1989). We include secondary electrons, the heat, 
and the Lyman alpha photons produced by primary electrons created via photoionization of 
HI, Hel and Hell by using results of a Montecarlo simulation (Shull & Van Steenberg 1985). 
No secondary electrons are considered when primary electrons are created via collisions 
with other atoms. This and the omission of collisional ionization of helium atoms are 
good assumptions since we consider only temperatures below 17,000 K. Unfortunately 
the complete set of results by Shull & Van Steenberg (1985) is given as function of the 
fractional ionization of H only for X-ray energies well above 0.1 keV. For lower energies 
they show the fractional energy deposited as heat and the number of secondary electrons 
ejected from hydrogen only for a discrete number of cases. We use a subroutine which 
interpolates their results and computes also the number of secondary electrons ejected 
from helium atoms, $i/ e , assuming that &He/&H is the same as at higher energies. A 
similar assumption is used for the fractional energy of the primary which goes into Lyman- 
a excitation for hydrogen (which will be relevant in computing subthermal effects). We 
introduce the following symbols: 

fiiTx), f2(Tx) are cooling efficiencies due to impact of neutral hydrogen (/i) and electrons 
(^2) with heavier elements (the contribution of each element has been weighted by its 
abundance with respect to H) . 

f^Tx) is the cooling efficiency for free- free emission and HI recombinations and f^Tx) 

is the cooling efficiency due to HI impact with free electrons. 

Xfj, XHe, and xu e 2 are the fractional ionization of HI, Hel and Hell respectively. 

hHe = 0.1, he = 1-9 x 10 -4 are the helium and carbon abundance by number with respect 
to hydrogen. 

ct2 is the HI recombination coefficient excluding captures to the n=l level (Spitzer 1978). 
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®He ®He2 are the Hel and Hell recombination coefficient (Spitzer 1978). 

«c is the total recombination coefficient for CI which we assume to be equal to that of HI 

and Ec = 2eV is the heat released for each carbon ionization. 

7# is the HI ionization rate due to collisions with free electrons. 

iai, $>ai,a2 and E a \ are respectively: the numbers of photoionizations per second per atom 
of type al, the number of secondary electrons created by interaction of type a2 atoms with 
primary electrons coming from photoionization of type al atoms, and the heat released 
for each photoionization of type al atoms. The symbol al can be H, He, or He2 to 
indicate HI, Hel and Hell atoms, the symbol a2 can be only H or He because secondary 
electrons from Hell atoms are not considered. We compute these quantities at z x / 2 for 
isotropic background radiations described in Section 2, attenuated by the overlaying and 
underlaying layers of gas. 

tin — Pj {kTx(l + XH + hHe + h>HeXHe + h>He%He2 + h>c)} is the volume density of H atoms. 
n e ~ uh(xh + hHeXHe + ^He^He2 + he) is the volume density of free electrons. 
The ionization equilibrium equations for H and He are then: 

n e { a 2 -1H) = ^H(l + $H,H + $H,He) H : hHeiHe 

V 1 — xh ' 1 — xh 

X (1 + <f> He ,H + $He,He) + h He £ He2 (l + $>He2,H + $He2,He) = C (8) 

1 — xh 

U e X H eOtHe = (1 - %He ~ X H e2)^He(^ + $He,He) + ~~T ~~~~ £, H H , H e 

+XHe^He2^He2,He (9) 



n e X He 2^He2 = X He2 ^He2 (10) 

The right hand side of equation (8) gives £, the total number of ionizations per second for 
a neutral H atom. If the slab is optically thin to the extragalactic background eq. (8) gives 
C the total number of ionizations per second for an unshielded neutral H atom. For Case 
A with a lower cut-off at 1.5 keV ( ~ 1.5 x 10~ 20 H ionizations s -1 , while (~3x 10~ 17 
H ionizations s _1 for Case C with E c = 0.1 keV, s = 2.4 and / ~ 5. The ionization rate 
C for Case B is about the same as for Case C with s = 2.4 and I xs 1.2. In paper II we 
shall extend the photon spectrum down to the Lyman edge, 13.6 eV, which results in much 
larger values for (. 
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The energy equation below describes the balance between the cooling rate and the heat- 
ing rate per H atom due to photoionizations by extragalactic background photons or to 
additional local heat sources (E nr ): 

ran {(1 - x H )fi(T K ) + [f2(T K ) + x H h{T K ) + (1 - XH )h{T K )] — } 
L riff ) 

= (1 - X H )E H £h + h He (l - X H e ~ X H e2)E H e^He 

+ h He x He2 E He 2^He2 + h c E c a c (T K )n e + E nr (11) 

The left hand side is the cooling per H atom minus the heat due to carbon ionization. The 
right hand side is the heating per H atom due to photoionization of hydrogen and helium 
plus the non radiative heating per atom, E nr . To simplify the solution of the above set of 
equations we set xu e 2 = at low temperatures while at high temperatures we omit the 
cooling function f\. Given iVffj, P , and the spectral index of the background radiation we 
solve the above equations by expressing Xff e , xn e 2 and P/£ as functions of xh using three 
of the four equations given above. For each value of T K then we need to solve numerically 
one single equation of type f(xn) = 0. 

-Case A: due to the high energies of the incoming photons, and to the low intensity of the 
radiation (which for this case is fixed) the gas is optically thin with negligible fractional 
ionization of Hell. Solutions of f(xn) = are found assuming constant values of <& a i,a2> 
E a \ and P. Since the degree of H ionization is low, P is determined only by P and N HI . 
We calculate <E> a i,a2 and E a \ (which are functions of xh) with an iterative procedure. Tk 
is determined by the value of P/£ if the heating rate due to H and He photoionization is 
stronger than that due to carbon photoionization. 

-Case B: for each value of Tk we guess a value of E nr and solve iteratively f(xn) = 
to find self consistent values of xh, &ai,a2, E a i, and P. The attenuated photon flux (or 
£) required at Zi/ 2 can be calculated. But since for this case the incident photon flux is 
fixed, we know what £ should be for a given xh- We adjust E nr accordingly and repeat 
the calculation until we get the correct value of £. 

-Case C: the value of xh, found by solving f(xn) = 0, determines P/£ and therefore I. 
The correct solution is found after a few iterations because the spectral distribution of the 
extragalactic radiation changes when the parameter I varies with Tk- 
For each case we then determine the spin temperature, as described at the beginning of 
this Section. 
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4. SUBTHERMAL EFFECTS IN THE 21-cm EMISSION LINE 



If the ionizing flux and the pressure P are moderate or small, the spin temperature Tg is 
depressed below the gas kinetic value Tk and this subthermal effect increases the 21-cm 
absorption (see Section 5). If flux and pressure were extremely small, Ts could be depressed 
so drastically as to approach Tr, in which case the observed emission quantity Nb would 
be much less than the actual column density Nhi- As mentioned in the introduction, "HI 
edges" are observed at Nb ~ Ni and the question arises whether HI could be "hidden" 
outside such edges. In other words, could there be regions where Nhi is only slightly 
smaller than Ni but Nb << Nhi ? We now give two arguments, one observational and 
one theoretical, to show that this is not the case. 

We selected a few objects with low 21-cm surface brightness (Nb < 2 x 10 19 cm~ 2 ) observed 
as part of outer disks by CS (4 in M33 and 1 in NGC 4631) or as HVC by CTS. These are 
good candidates for a phenomenon of N HI » N B - In Figure 2 we show the lower limits 
to Ts calculated from the absence of 21-cm absorption (we assume a line width of 4 km/s 
except for one case in M33 where Ts > 300K and we use 8 km/s). The curves give the 
Nhi — Ts relation (see eq. (4)) for each selected object and the thick part of the curves 
displays the allowed values of Nhi- Since they are all on the vertical portion of the curve, 
where Ts >> Tr and Nhi — Nb, we can exclude that strong subthermal effects hold for 
these cases. 

For the lowest extragalactic flux and pressure, namely Case A with E c = 1.5 or 0.5 keV and 
Po/k = 0.1 cm -3 K, we have calculated Tg. These small values of flux and pressure should 
give the largest difference between N B and Nhi- Figure 3 shows the N B —N H i relation and 
we can see that no drastic depressions in the brightness are expected for Nb > 10 18 cm -2 
(appreciably smaller than Ni). Therefore this proves again that the surface brightness is 
still a good indicator of the neutral hydrogen column density where HI edges occur. Our 
conclusion that HI edges are not due to subthermal effects holds even more strongly if the 
background spectrum extends down to lower energies. Table 1 gives Nhi/Nb, Tk and Ts 
for Case A with E c = 1.5 and 0.5 keV and for a few low values of Pq. Notice that the 
gas is always in the cold phase and as we increase the pressure Ts approaches Tk, Tk 
decreases until reaches 12.5 K at very large Pq. This is the lowest possible temperature 
for our assumed metal abundance and it is reached when heating is provided by carbon 
ionization. For larger ionizing fluxes (e.g. Case C) we have larger values of Tk and of Ts 
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with a consequent decrease of subthermal effects in emission. Values of Ts in Table 1 are 
already too large to allow substantial departures of Nb from Njji, but they are smaller 
than those observed by CS and CTS (some of which are shown in Fig. 2). This implies 
extra heat sources in outer regions (Case B) or larger ionizing fluxes than we used for Case 
A (case C). 



5. THE SPIN TEMPERATURE - HEAT RELATION 



The flux assumed in Case A is unrealistically small and implies spin temperatures well 
below those observed. We confine therefore our attention to Case B and Case C which 
involve stronger energy inputs to the gas than Case A. In Case B the ionizing flux is fixed, 
the non radiative energy input rate E nr varies, and Tg (as well as T K ) are calculated. In 
Case C, E nr = and the intensity of the ionizing flux for E > 0.1 keV, /, varies. The 
results we show are useful for interpreting present and future HI absorption data in terms 
of heat inputs to the gas from local sources or, in the absence of these, from a cosmic 
background flux. Using the HI data we have available at the moment, we show that the 
heat input required from local sources, E nr , is small and compatible with the absence 
of star formation in outer regions. If there are no local sources (Case C) a comparison 
between the observed and predicted spin temperature for a given HI column density sets 
interesting constraints on the intensity of the cosmic background radiation. 
We know that for given P and Tk a medium with a lower ionization fraction in the cold 
phase has a lower cooling rate per atom (Spitzer 1978). This means that for a given T K 
Case B requires a heating rate per atom (radiative plus non radiative) lower than Case C. 
However, because subthermal effects in Case B are stronger than in Case C, this inequality 
might not hold for a given T$. Our results show that even when the external pressure is 
low and subthermal effects are strong, there is relatively little difference between the heat 
required in Case B and in Case C for a given Tg. Figure 4 plots E nr versus Ts for Case 
B in the stable cold and warm phases (solutions requiring E nr < have been discarded). 
For Case C, Figure 5 plots the parameter / in eq. (3) versus Ts for s = 2.4. Notice that 
for very small values of I the extragalactic radiation field is unimportant in the energy 
equation, Ts is constant and determined by the balance between the carbon heating and 
the line cooling. In Figure 6 we plot Tk versus Ts for all the cases shown in Figure 4 and 
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Figure 5. Notice that the warm phase can starts with lower spin temperatures than the 
final part of the cold phase. 

Figure 4 and 5 show E nr or / versus T$ for three different values of N HI and two values of 
the additional pressure Po- As we increase the parameter P , and therefore the density, the 
heat input (E nr or I) required to reach a certain Tk increases but subthermal effects get 
smaller and Ts closer to Tk- The result for a given Ts is that E nr or / increases rapidly 
with P if Ts is low. For high Ts the heat input required is almost independent of P and, 
if the density is not big enough that Ts stll differs from T K , it can increase with P (see 
for example Fig. 5 for Ts ~ 1000K). 

For Case C the existence of the cold and/or warm phase of HI depends mainly on the 
ratio of P (from eq. (7)) to the intensity /: for P/I < 3k there is no cold phase at all, for 
slightly larger values both phases coexist, and for P/I > 10k there is no warm phase. If 
the contribution to / from known quasars alone gives / ~ 5, the outermost regions of the 
HI disk, with column densities Nhi ~3x 10 19 cm -2 , are likely to be all in the warm phase 
when Z1/2 ~ 1 kpc. This should be true even in the absence of E nr or other additional 
ionizing fluxes (a warm phase with smaller zi/ 2 requires stronger background fluxes). This 
value of Z1/2 is what we estimate from Merrifield (1992) for the outermost disk of the 
Milky Way. Therefore even if the outermost part of galactic HI disks are warm there is 
no strong evidence for the heat input in outer regions to be much larger than that given 
by Madau's flux. Figures 4,5 and 6 show a curious feature for Case B with low heat input 
(E nr ^ 10 -15 eV s _1 ) and for Case C with / ?S 10: the hydrogen is all in the warm phase 
with Tk ~ 10 4 K but the spin temperature is surprisingly low, Ts 500K. Conversely, if 
observations should give Ts ^ 1000K much larger ionizing fluxes and/or larger E nr would 
be required. 

We illustrate the use of Figures 4, 5, and 6 by applying them to the analysis of the 
observational data for the strongest background source behind the outer disk of M33 (the 
curve farthest to the right in Fig. 3). Here we have Nhi — Nb — 2 x 10 19 cm -2 along 
the line of sight, and the absence of absorption gives Ts > 250 K (Corbelli & Schneider 
1990). The HI surface density perpendicular to the plane will be slightly smaller than the 
observed value due to the inclination of the outer disk respect to our line of sight (which 
is less than that of the inner disk but not completely negligible, see Corbelli et al. 1989). 
Fig. 4 (b) and 5 (6) show that we have possible solutions for both Case B and C with small 
heat input if Pq is reasonably small (Po/k ^ 10, say). The minimum heat compatible with 
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Ts ~ 250 K for small pressures is similar for Case B and C and gives solutions in the 
warm phase with E nr ~5x 10 -16 eV sec -1 for Case B and I ~ 6 for Case C with s = 2.4 
and E c = 0.1 keV. This value of / can be achieved by quasars alone (Madau 1992) and 
this value of E nr can be supplied by a mild outer galactic fountain (Corbelli & Salpeter 
1988). For instance, one requires only 1% of the total galactic supernova energy output 
rate to flow with energy of ~ 10 5 eV cm -2 s -1 over a disk of 30 kpc radius. If half of 
this energy is used for heating the disk gas layer below, then Ts ~ 500 K for Nhi ~ 2 or 
3 x 10 19 cm -2 . Slightly smaller spin temperatures are predicted for a gas with lower N HI . 
To summarize: the present limit of T$ > 250K is still compatible with present estimates of 
the quasar flux and/or a very mild outer galactic fountain, which are sufficient to keep a 
gas with N HI < 3 x 10 19 cm -2 in the warm phase. The spin temperatures are surprisingly 
low, compared to Tk ~ 10 4 K (unless the energy input is extremely large). 
Another reliable absorption measurement of 21-cm line outside the optical disk is in the 
spectrum of the quasar 3C232, at the velocity of the spiral galaxy NGC3067, close to 
the quasar on the plane of the sky. The width of this line (Rubin et al. 1982) implies 
that the absorbing material is HI in its cold phase with Ts < 300 K and with N^f\ d < 
5.8 x 10 19 cm -2 along the line of sight as derived from the strength of the absorption line. 
Emission measurements (Carilli & van Gorkom 1992) give a total observed HI column 
density of (8 ± 4) x 10 19 cm -2 ; assuming a modest inclination correction factor, it should 
be Nhi ~ 5 x 10 19 cm -2 , so that Figures 4(c) and 5(c) apply. We can be in the narrow 
range of conditions where the warm and the cold HI phase coexist. If Zi/ 2 is of order 
1 kpc, as observed in the Milky Way, then Po/k ~ 30; in the vicinity of this pressure 
Fig. 5(c) gives a flux with / ~ 3(P /I0k) - 5 and Figure 4(c) gives for Case B a heating 
rate E nr ~ 5 x 10~ 16 (Po/10/c) 3 eV s _1 . These are only order of magnitude relation, 
but they are compatible with the estimates given above for M33 and with likely quasar 
fluxes. Figure 4(c) and 5(c) also make predictions about the warm phase: Ts should be 
~ 300 — 500 K for Case B and only slightly larger for Case C. More than half of the HI 
should be in this phase and one might be able to observe a second absorption component 
with the peculiar combination of a low spin temperature and a large width, u ~ 20 km 
s _1 , corresponding to the larger Tk- 
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6. DARK MATTER DECAY AND HI TEMPERATURES IN OUTER REGIONS 



The dark matter decay theory recently developed by Sciama (Sciama 1990, 1991), predicts 
a universe populated by heavy neutrinos which decay radiatively producing photons with 
energy E 1 = 13.6 + e with e ~ 1 eV and a production rate proportional to the density of 
dark matter. If the local density of the gas is sufficiently high, these photons are produced 
and absorbed locally and they are the most important photons for the ionization of HI. 
In the outermost regions of galactic disks the local conditions are very different and the 
monochromatic photons, both extragalactic and local, can cause a sharp HI edge, even 
in the absence of quasar photons. This is discussed in Paper II. Here we again consider 
only the mostly neutral portion of the galactic disk, just inside the edge. We assume first 
a negligible background flux from quasars and no heat source besides the local neutrino 
decay photons. Due to the larger number of photons with energy close to 13.6 eV, the 
Lyman-a pumping will be very effective to bring T$ close to Tk- As source function for 
these photons we use the expression given by Sciama (1990) for our Galaxy: 

5 x 10" 16 o i 

^ 3x(l + W Cm S <12) 
The equation of state relates density and pressure; the ionization recombination equilibrium 
and the energy equation in the absence of any external flux are: 

PXH , r (13) 




T K (l + x H + h He ) 



T K (1 + x H + h He ) LV a 2 T K (l + x H + h He ) J 

Combining these two equations we have a direct relation between P and Tk, which depends 
on the square root of the source function 0, and on e through xh- For R/Rq = 3 Figure 
7(a) shows this relation and Figure 7(6) the corresponding fractional ionization xh- Notice 
that the gas with low HI column density can survive the external flux in the neutral state 
only if its pressure is higher than 4 cm -3 K; there is also a range of pressure which allows 
two different equilibrium temperatures. Absorbing the external flux, F ext due to neutrino 
decay in the halo and in the rest of the Universe, requires a column density larger than 
N m i n = 2F ext TK I '(P0L2), which is plotted as the thick curve in Fig. 7(c). The thin curve in 
the same Figure shows the maximum column density compatible with a given Tk (obtained 
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by setting Pq = 0). Therefore in the absence of any other heat source Tk ^ 500 K and the 
region between the two curves indicates all the possible Nhi — Tk pairs. 
If there are extra heat sources in the medium, Tk can be higher than the values shown in 
Figure 7 and the energy equation will be eq. (11) with the helium terms neglected and the 
heat from a quasar background replaced by the heat from the local monochromatic flux. 
The non radiative heat input which is needed is: 

P(l-x H ) , , r, , , , wli ,T K (l + x H ) 



Enr — m , x /l + 

T K (1 + x H ) 



h + x H h + (1 - x H )U - - ^ K p (15) 



where xh = a/'/W ^/ {PjT K — y l 4> v /ot2)- E nr depends on Nhi only through the pressure 
term and therefore we can plot E nr versus T K for several values of P without using P 
and N H i explicitly. Figure 8(a) shows the E nr - T K relation for P/k = 10, 100, 1000 
cm -3 K. The extent of the curve is limited by the conditions E nr > 0, < xh < 1 and 
Nmin < N H i < N max (N min and N max being defined as in Figure 7(c)). For the same 
values of P, Figure 8(6) shows N min (thick curves) and N max (thin curves) as functions of 
Tk- Suppose we have gas with Nhi — 2 x 10 19 cm -2 and T$ > 500K; then the minimum 
required extra heat will be E nr ~ 2xl0 -15 eV s _1 which is comparable to the E nr required 
for the same column density by Case B, where there are no local monochromatic photons 
but a background at higher energies. For this case the minimum value of P/k needed for 
HI survival against the neutrino's external flux is ~ 10 cm -3 K. Higher column densities or 
higher pressures require slightly higher values of E nr . These E nr values are in reality upper 
limits since even a modest quasar background, which we have neglected, will contribute 
somewhat to the heat. 

We return now to the two examples of spin temperature measurements. For Nhi ~ 2 x 10 19 
cm -2 as in M33, we only have the observational limit T$ > 250K. Figure 7(c) shows 
that this is compatible with neutrino decay theory even in the absence of extra heat 
sources. However for a surface density N H i ~ 5 x 10 19 cm -2 (as it might be in the case 
of NGC3067/3C232) Figure 7(c) gives T s < 30K (T K ^ T s ) whereas the observations 
showed a higher spin temperature and probably some warm HI. Non-ionizing heat sources 
are therefore required, but E nr need not be larger than for Case B. 
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7. SUMMARY AND DISCUSSION 



In this paper we have discussed HI emission and absorption studies at 21-cm wavelengths 
which can be used to constrain the spectrum of the extragalactic ionizing radiation below 
1 keV, and more generally the rate of heat input outside the optical disk of galaxies. 
Taking into account self-gravity and considering only unreasonably small lower limit to the 
extragalactic flux (Case A) we have shown that the number of ionizations and collisions are 
sufficient to bring the spin temperature well above the background temperature Tr even 
in the absence of strong compression (due to coronal gas or dark matter) . This means that 
the actual column density, Njji, does not differ appreciably from the value Nb, inferred 
from the 21-cm brightness temperature at the periphery of galaxies. HI edges observed in 
outer disks around Ni ~ 2 x 10 19 cm -2 (Corbelli et al. 1989; van Gorkom 1991), where 
Nb drops very rapidly below the minimum detectable value, or the small quantity of HVC 
gas with N HI < 5 x 10 18 cm -2 , correspond then to real cut-offs in the neutral phase of 
the hydrogen distribution and not to a sudden decrease of the HI spin temperature. These 
edges are most likely due to ionization by the background radiation and occur when the 
total gas column density drops below a critical value. This phenomenon will be analyzed 
in Paper II. 

We have carried out model calculations for both the gas kinetic temperature, Tk, and the 
HI spin temperature, T$, for various values of gas density and extragalactic background 
fluxes. We consider a fairly narrow range of hydrogen column densities appropriate to 
outermost HI disks of spiral galaxies, Nhi ~ (0.5 to 5) xlO 19 cm -2 , which are most 
sensitive to photon energies ~ 0.1 keV. We consider fluxes from below to above the current 
range of estimates of quasar backgrounds at zero redshifts. As in the classical "two-phase 
model" we find for these backgrounds T K ~100K for the cold phase and T K ~ 10 4 K for the 
warm phase. Deviations from thermal equilibrium are still strong and T$ « T K although 
Ts » 3K always. We have carried out calculations both with and without additional 
non-ionizing sources, E nr (to model mild galactic fountains, hydromagnetic waves, etc., 
traveling from an inner disk to an outer disk). The relationships between E nr , N HIl Ts, 
Tk and the intensity I of the background flux are given in Figures 4, 5, and 6. 

The models depend not only on the intensity / (for E nr = 0), but also on an additional term 
Pq, with which we estimate the effective compression of the gas (dark matter compression 
versus expansion due to some internal non thermal pressure) . If we assume that outermost 
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HI disks with Nhi ~(2 to 5)xl0 19 cm" 2 have a vertical scale height Zi/ 2 ~ 1 kpc, as 
found for our Galaxy (Merrifield 1992), we obtain Pq ~ (10 to 30) cm" 3 K. Our models 
predict that the HI gas in low column density regions should be in the warm phase but the 
gas in regions with column density of order 5xl0 19 (which have a slightly higher volume 
density) should be a mixture of cold and warm phases. We analyzed data from an outer 
region of M33 with Nhi ~ 2 x 10 19 cm -2 where only an upper limit to the absorption 
was found (Corbelli & Schneider 1990) and one region in NGC3067 with N HI ~ 5 x 10 19 
where some absorption was detected (Carilli & van Gorkom 1992). The non detection of 
21-cm absorption lines through the outer disk of M33, which is a relatively undisturbed 
HI disk, points out that a significant fraction of the HI gas is warm despite the absence of 
star-forming regions (in the inner, star-forming disk the fraction of cool HI for M33 is only 
slightly smaller than for the Milky Way; Dickey & Brinks 1993). This requires a cosmic 
background stronger than given by the extrapolation of the hard X-ray spectrum down to 
0.1 keV. The absorption/emission study for NGC3067 suggests that both the warm and 
the cold HI phases are present at the slightly larger column density there. 
We found that both observations are compatible with an intensity / of extragalactic pho- 
tons at energies ~ 0.1 keV roughly as estimated by Madau (1992) for the present day 
quasar background. If this extragalactic ionizing flux should turn out to be smaller, one 
can still explain the observed temperatures by adding a modest amount of non ionizing 
heat input E nr . Similar heat inputs are required if one replaces the external ionizing flux 
by monochromatic UV photons produced locally by neutrino decay (Sciama 1990). Only 
if future observations will show that HI spin temperatures at the periphery of a galaxy 
are as large as 1000 K, would strong local heat sources be required in addition to a QSO 
ionizing background. 

The warm HI phase of the ISM is a dominant component in outermost regions and we can 
make the following prediction: besides possible narrow HI absorption lines from any cold 
phase, the warm HI should have a spin temperature still very different from the kinetic 
one. Values of Ts ~ 500K should be common unless there are heat inputs, either from a 
background or from local sources, much stronger than what is now believed. This warm 
phase should produce an absorption line component with quite large width (a; ~ 20 km 
s _1 ), which is difficult to detect but would provide an interesting diagnostic because of the 
Ts — iv mismatch. 
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FIGURE CAPTIONS 



Figure 1. Some of the X-ray spectra used. At energies above 1.5 keV all the spectra have 
s = 1.4 / = 7.7, and the thick line shows this spectrum extended down to 0.1 keV as 
it has been used for Case B. The same spectrum has been used for Case A but with 
lower cut off energies at 1.5 or 0.5 keV. All the other lines are examples of spectra 
used for Case C and are labelled according to the corresponding values of / and s 
(see eq. (3)). 

Figure 2. The Ts — Nhi relation for all the positive detections of HI with Nb < 2 x 10 19 
cm -2 in CS (circles) and Nb < 10 19 cm -2 in CST (triangles). The thick parts of 
the curves show the possible pairings of Nhi and Ts for each data point according 
to the upper limits obtained from 21-cm absorption data. 

Figure 3. Subthermal effects in emission: the column density derived directly from the 
brightness temperature, Nb, versus the real Nhi of the gas. We have used Po = 0.1 
cm~ 3 K, and a low external flux, namely Case A with E c = 1.5 keV (dotted line) or 
E c = 0.5 keV (dashed line). The solid line is the line of unity slope in the Nb-Nhi 
plane. 

Figure 4. Heat per atom required for equilibrium in units of eV s -1 for Case B (E c = 0.1 
keV, s = 1.4, / = 7.7). It is shown as function of Ts in the stable cold and warm 
phases for two values of Pq- The filled triangles indicate the beginning of the two 
phases when Po/k=l cm -3 K (solid lines) and the open triangles when Po/k = 100 
cm -3 K (dashed lines). N H i = 5 x 10 18 cm -2 in (a), 2 x 10 19 cm -2 in (6), and 
5 x 10 19 cm -2 in (c). 

Figure 5. The parameter / of the spectrum in eq. (3) as a function of Ts for Case C 
when s = 2.4. The filled triangles indicate the beginning of the cold and warm phase 
when P /k=l cm -3 K (solid lines) and the open triangles when Po/k = 100 cm -3 K 
(dashed lines). Nhi is 5 x 10 18 cm -2 in (a) is 2 x 10 19 cm -2 in (b) and is 5 x 10 19 
cm -2 in (c) 

Figure 6. The kinetic temperature Tk-Ts relation for all cases shown in Fig. 4 and in 
Fig. 5. (a), (6), (c) are for Case B with Nhi = 5 x 10 18 , 2 x 10 19 , and 5 x 10 19 cm" 2 
respectively. For the same values of Nhi (e), (/), (d) show the Tk-Ts relation for 
Case C with s = 2.4 and E c = 0.1 keV. Lower curves in the panels are for the cold 
phases while curves at the top of the panels are for the warm phases. 



23 



Figure 7. In (a) and (b) we show the pressure and the hydrogen fractional ionization for a 
gas pervasive of monochromatic photons of 14.6 eV from decaying neutrinos and no 
other heat sources. The thin and thick curves in (c) show the maximum and minimum 
HI column density, respectively, which are compatible with the corresponding value 
of Tk (see text for details). 

Figure 8. In (a) the extra heat per atom required by a gas with P = 10 (continuous line), 
P = 100 (small dashed line) and P = 1000 (large dashed line) cm -3 K, is shown as 
function of the kinetic temperature. In (b) thin and thick curves are the maximum 
and minimum HI column density respectively, which are compatible with a given 
value of T K (see text for details). 
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